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Abstract 

In this paper, a mathematical model on the connecting hose of a one-stage Oxford split-Stirling cryocooler has been made to 
simulate its dynamic performance. A set of governing equations has been discretized by method of characteristics and dissolved by a 
numerical method. The present paper emphasizes the validation of this model against experimental results of pressure phase lag and 
amplitude attenuation through connecting hose, and cooler performance vs. connecting hose in different dimensions. The com¬ 
parison shows acceptable agreement between the simulation and experimental results. 
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1. Introduction 

The Stirling cycle cooler comprises a compressor and 
a cold finger, which are connected by a connecting hose. 
This cooler operates on a closed regenerative cycle by 
fixed gas mass being cycled repeatedly through the sim¬ 
ilar thermodynamic sequence; a moving displacer/re¬ 
generator in the cold finger and a compression piston in 
the compressor are driving by a linear motor respectively 
with a fixed phase difference between them; the thermal 
regenerator is the key, central component. Comprising 
nothing more than a porous, conducting matrix with a 
high ratio of wetted area to volume, this passive, de¬ 
ceptively simple component stores and releases heat in¬ 
ternally which would otherwise be exchanged with the 
surroundings at temperatures intermediate to the cycle 
limits of T c and 7 e. To this extent, the closed regenerative 
cycle fulfills the basic condition for thermodynamic re¬ 
versibility. 

As we well have known, miniature cryogenic cryo- 
coolers operating reversed Stirling cycle are capable of 
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operating over a range of temperature differentials, has 
been produced commercially in a wide range of heat lift 
capacities, and are used for cooling the detectors of 
thermal imaging systems. The type is a candidate for 
duty in satellites—a singularly demanding environment, 
since whatever the cooling power achieved, there is 
pressure for it to be greater; whatever the electrical 
power consumption, there is demand for it to be less. 
The Oxford split-Stirling-cycle cryocoolers, which are 
having been developed since its invention by the re¬ 
searchers in Oxford University at the end of 1970’s [1], 
and many special components or technologies such as a 
linear compressor, no-contact seals, flexure springs, 
controlling contaminants technologies and etc., have 
been introduced during the manufacturing and assem¬ 
bling processes, and is being one quite popular kind of 
miniature variants in many fields, space applications in 
particular, tend to operate between 60 K and ambient, 
lifting a fraction of a watt for an electrical input of a few 
watt [2,3]. 

An Oxford split-Stirling-cycle cryocooler is com¬ 
monly made up of a linear compressor and an expander, 
which are connected by a flexible hose. The regenerator 
always contains a long stack of phosphor bronze gauze 
discs. The compression piston and the displacer move at 
a high frequency (30-60 Hz), so the cycle of work, heat 
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Nomenclature 


A heat transfer area 

A x cross area of connecting hose 

A w wall area 

a local sound velocity; constant in Eq. (5) 

b constant in Eq. (5) 

c p gas specific heat at given pressure 

c v gas specific heat at given volume 

D diameter 

D eq equivalent diameter 

D h hydraulic diameter 

/ friction factor, fining factor or frequency 

h heat transfer coefficient 

k thermal conductivity 

m gas mass 

m mass flow rate 

p pressure 

Q heat flux rate 

Q+,Q- discretization factor 
R gas constant 

Re Reynold’s number 

S stroke 

5 entropy 

T temperature or period 

7j , T discretization factor 

u gas flow velocity 

V volume 

x coordinate of flow direction 

Greek letters 

a medium parameter 


friction parameter 
5 medium parameter 

e convergent criteria factor 

<I> differential phase angle 

A characteristic factor 

/< dynamic viscidity or modified empirical co¬ 
efficient 

9 representative parameter 

p density 

a convergent constant 

t time 

y ratio of specific heat 

energy parameter or modified resistance co¬ 
efficient 

Subscripts 

C compression chamber or compressor 

D displacer 

g gas 

h hydraulic 

H hot chamber 

max maximum 

new new value 

old old value 

P large section of the assumed nozzle 

T small section of the assumed nozzle 

t any section of straight tube or hose 

W wall 


exchange and fluid flow variations is a complex dynamic 
and thermodynamic interaction. In the hands of Fin- 
kelstein [4], Organ [5], Huang [6], Yuan [7,8], Sun [9,10], 
de Monte [11] and etc. the performance simulation and 
experimental verification have been taken to a high art 
form. On the other hand, because of most of the relevant 
phenomena being, in reality, non-linear, and while lin¬ 
earisation of an isolated, non-linear phenomenon may 
be justified, no-one knows for certain how much extent 
of performance irreversibility of reality could be intro¬ 
duced by the connecting hose. The connecting hose 
could give rise to phase lag and amplitude attenuation of 
gas pressure wave and mass flow rate wave, and what’s 
more, this could change the phase relations between the 
mass flow and pressure oscillation in different chambers, 
which are one of the important factors that influences 
the performance of regenerative cryocoolers. Unfortu¬ 
nately, about the connecting hose’s operating charac¬ 
teristics and its effect on the cooling performance of the 
Oxford split-Stirling-cycle cryocooler, there are so few 


reports or publishing results at present. All these make it 
very difficult for the researchers to grasp the genuine 
performance characteristics, and rely on a lot of ex¬ 
periments or experimental data during the thermal 
design process, the connecting hose’s operating charac¬ 
teristics and its effect on the cooling performance of the 
Oxford split-Stirling-cycle cryocooler have been put 
forward by us and studied extensively. 

To the best of the author’s knowledge, a mathemat¬ 
ical model for the connecting hose of an 80 K Oxford 
split-Stirling-cycle cryocooler has been described in this 
paper. A set of differential equations based on the con¬ 
servation of mass, energy, and state has been discretized 
by method of characteristics and dissolved by a nu¬ 
merical method. Furthermore, the cooler performance 
validation for the cooler with connecting hoses in dif¬ 
ferent dimensions against experimental results for an 80 
K Oxford split-Stirling-cycle cryocooler have been per¬ 
formed and acceptable agreement was found between 
the two. 
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2. Gas process model for connecting hose 

2.1. Physical model and governing equations 

The physical model for an Oxford split-Stirling 
cryocooler is shown in Fig. 1. In order to set up a 
comparatively practical mathematic model, the basic 
assumptions are made as follows: 

1. one-dimensional flow and heat transfer; 

2. neglecting the length of the trumpet inlet or outlet at 
both ends of the connecting hose; 

3. ideal gas and constant properties; 

4. negligible axial conduction along the connecting 
hose; 

5. negligible clearance leakage between the cylinder and 
piston, or the cylinder and displacer; 

6. constant wall temperature of the hot and cold cham¬ 
bers. 


The following describe dynamic model and governing 
equations for working gas in the connecting hose [12], 
Continuity equation: 


OP 

0T 


0w dp apu 

p- + Ul f + -!—=0 

0X 0X X 


( 1 ) 


where a = (d lm4 v /dlnx) = x(dlm4 v /dx). Note that the 
mass flow velocity u is defined as positive if the gas flows 
form the compression chamber to the hot one and 
negative for the opposite direction. And so it is with the 
mass flow rate, m = pu(nD\)/A. 





Fig. 1. (A) Physical model; (B) regarding inlet or outlet as nozzle. (1) 
Compression piston, (2) compression chamber, (3) connecting hose, (4) 
hot chamber, (5) displacer/regenerator. 


Momentum equation: 


0M 

i> Ft + P „ 


du 

dx 


dp 

0X 


P 


( 2 ) 


where = — (8f /D h )(pu 2 /2)(u/\u\), f is the fining fric¬ 
tion factor, D h is the hydraulic diameter of the con¬ 
necting tube. 

Energy equation: 


dp dp , 
2-+«2--a 2 
0T dx 


dp dp'. 

£ + “£ )=* 


( 3 ) 


where i p = (y — 1 )(puaQ x - up), y = c p /c v , Q x is the heat 
flux rate at the position x. 

Equation of state: 

P = PRT ■ (4) 


2.2. Instantaneous heat transfer coefficient and the friction 
factor 

Values are required of the coefficient of convective 
heat transfer, h, for use in the energy equations (Eq. (3)). 
The task is best done in terms of Stanton number, and 
corrected in turn with Reynolds number and Prandtl 
number for conditions in the variable-mass or variable- 
volume system. The only source known to the authors is 
An Hand's correlations [13] for conditions in the cylinder 
of the reciprocating internal combustion engine. 
Annand determined values of a and b for the following 
expression for heat flux rate: 

Q = - Tg) + eo 0 (l* - T g 4 ) J (5) 

in which the Reynolds number, Re, is defined as 
Re = 2fSD/p, where S and D is the stroke and diameter 
of the compression piston, respectively; / is the piston 
motion frequency. As for the connecting hose, we con¬ 
vert an equivalent D eq = LD h /S, where L and £> h are the 
length and hydraulic diameter of the connecting hose, 
respectively. 

Neglecting the radiation component as being insig¬ 
nificant at temperatures involved, and defining h as the 
ratio of the heat rate, Q, to the product of wetted area, 
Aw, with temperature difference, 7 W - T g , leads to 

h=a^Re h (10 3 <Re< 10 5 ) (6) 

As for the fining friction factor f, the following empir¬ 
ical formula for oscillating gas flow in long tube from 
Roach and Bell [14] was adopted. 

/ = 0.1556[p t ( Mt ) ma ^ th M]- a201 [( Mt ) max Mh] (7) 

in which (« t ) max an( l ^th are the maximum velocity of the 
oscillating gas flows, and the hydraulic diameter of the 
connecting hose respectively. 

It is worth noting that / is considered to be a constant 
during the modeling, but should be adjusted by nu¬ 
merical iteration during the computation to give a 
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correct value for the responding maximum gas velocity 
(wt)max- For a given forward node i, p t and m, is known, 
Pi+i and u i+ 1 of backward node i + 1 can be determined 
by using m, as mean velocity of this segment and / from 
Eq. (7). Then we use the average of m, and u i+l as mean 
velocity and / from Eq. (7) to determine (p/+i) new If the 
two values of p i+ \ is not satisfying to Eq. (8-1), then a 
new / can be determined from p- t and p i+ \. Comparing 
new / with / from Eq. (7), we could determine the ad¬ 
justing value and symbol till the value of / satisfies the 
criteria of the following two equations: 

l(A+i) new (P-+i)oidl < £ ^ o (8-1) 

l(A+l)newl 


/ 


pc&Vc - f 


Ph d Fh — Wk — Qh 


Qx 


<b 


d P 


dp 

0X 


i ^P j 

dx + f- dt 

0 T 


( 11 ) 


Now we can set out a set of six, simultaneous equa¬ 
tions for the property gradients in systematic form, then 
could gain the solution for any partial coefficient ex¬ 
pressed in determinant form. The mathematical condi¬ 
tion that the gradient should be discontinuous is that the 
determinant expressing the gradient should be indeter¬ 
minate—i.e., should be 0/0. Therefore, we could get the 
following expressions for the local gradients of the 
Mach lines and the characteristics line (or path line) as 
follows respectively: 


r) = 

d x J ± u± a 


( 12 ) 


=> 0 ( 8 - 2 ) 

where ffjf represents the friction work in one cycle of gas 
friction in the connecting hose; Q\\ and Q t represent the 
heat loss in one cycle of gas in hot chamber and the 
connecting hose to the environment, respectively; s and 
d can be adjusted according the solution accuracy and 
computational speed. 


, ^ d P 

du = a - 

P 


P a dT 
p T dx 


uacoc 


dt 


(13) 


The compatible equations for path lines and Mach lines 
are given, respectively, as follows: 


2 I fl 

dp — a dp = — dx 
u 


(14) 


2.3. Process irreversibility 

For the irreversibility or entropy generation caused 
gas flow through the connecting hose, we utilize a 
standard thermodynamic expression for change in en¬ 
tropy corresponding to changes T and p, which may be 
written [15] 

T ds = c r dT + pd(\ / p) (9) 


)d t± (15) 

Gas physical or Mach line for connecting hose is shown 
in Fig. 2. Eqs. (14) and (15) can be converted the fol¬ 
lowing expressions by using finite difference format re¬ 
spectively [16] 

P4 + Q+u 4 = T+ (16-1) 


d p± ± padu± = ( — c P ua _|_ a p _|_ ^ 


For a process involving heat transfer and friction over a 
time interval dr the change in entropy, d,v. is 


ds = 




( 10 ) 


p 4 - Q-U 4 = T_ (16-2) 

where Q .Q ■ T .. T are parameters adopted during 
calculating process. 

The following section will discussed the dealing 
methods or measures to the boundary nodes. 


3. Dealing by method of characteristics and boundary 
conditions 

3.1. Discretisation of governing equation 

We proceed by regarding Eqs. (l)-(3) as equations in 
the gradients, ow/ox, du/dr, dp/dx, dp/d t, dp/dx, du/dr, 
but so far only three differential equations. If, however, 
the assumption is made that the fluid properties them¬ 
selves, u,p , p are continuous over the x — r region to be 
considered, then use can be made of the three further 
differential expressions which define this continuity. The 
expression for pressure typifies those for u and for p. 
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3.2. Dealing with the boundary conditions 

3.2.1. Gas pressure and volume variations in compression 
chamber 

For the one-stage Oxford split-Stirling cryocooler, 
the pressure wave is provided with a linear compressor 
coupled to the hot chamber through a connecting hose, 
as shown in Fig. 1. This design can make the coldtip not 
be affected greatly by vibration from the compressor. 

Since the compression piston motion is very fast and 
the gas mixing in the compression chamber is severe, the 
gas temperature and pressure inside the compression 
chamber are assumed to uniform. The compressor cyl¬ 
inder is assumed to be adiabatic in the analysis, since 
each of the compression and expansion processes occurs 
in such a short period of time that little heat exchange 
between the gas and the cylinder wall could be effectu¬ 
ated. The gas adiabatically compressed in the cylinder is 
assumed to be cooled to room temperature by the cyl¬ 
inder lid, compressor base frame and the connecting 
hose, so that the temperature of the gas entering the hot 
chamber is always equal to room temperature. 

Applying the first law of thermodynamics to the 
control volume drawn around the volume swept by the 
compression piston in the cylinder, the following equa¬ 
tions can be obtained for ideal gas 

d V c = — (rTz d m c - — dpc^l (dm c >0) (17) 

pc \ y J 

dV c = Vc(—~y & ^-\ (dm c <0) (18) 

V me Pc J 

If m c is defined as positive when the gas is leaving the 
cylinder, the conservation of mass gives 

d m c = —m c dr (19) 

For the Oxford split-Stirling-cycle cryocooler, the 
compression piston is driven by a linear motor. The 
equation of volume variation for the compression 
chamber can be written as follows: 

V c =Vc + AV C sin(27t/r) (20) 

where V c and AV C are, respectively, the average volume, 
and volume variation amplitude of the gas in the control 
volume; and / is the frequency for the moving com¬ 
pression piston. In practice, volume variations for the 
compression chamber and cold chamber are not exactly 
sinusoidal wave; so their equations should be rectified 
according to their measured results. 

State equation for the gas in the compression volume: 

p c V c = m c RT c (21) 

From the above Eqs. (17)—(21), the variation with time 
of p c and m c . may be determined by numerical calcu¬ 
lation. 


3.2.2. Gas volume variations in hot chamber 

The displacer is driven by a linear motor, and the 
motion phase difference maintains constant ( I> with the 
compression piston motion. The equations of volume 
variation for the hot chamber can be written as follow: 

V H = V H + AV H sm{2Tifx + n + <P) (22) 

Being similar to Eq. (20), Eq. (22) should be rectified 
according to experimental results [2]. 


3.2.3. Gas characteristics in the nozzle 

The gas characteristics in the nozzle can be described 
by the following governing equations: 

Continuity equation: 


OmpPp = u T p T 

Energy equation: 


2 2 , y ' 

a„ — «p H- 


1 


Up - Qt 


y + 1 


a p = 


Local sound velocity equation: 

jypp a _ jyp t 

Pp V Pt 

Isentropic process equation: 

Pr = (PrY = ( a rV y/ ^ 

Pp V Pp / V «p / 

Momentum equation: 

OppMp - p T Uj = pr - Pp 


(23) 

(24) 

(25) 

(26) 

(27) 


No matter which direction the mass flow is in, one 
nozzle acts as inlet and the other acts as outlet. The 
below equations always could come into existence. 

Pp = Pc or pp = p H . (28) 


4. Calculation method and procedure 

From the above description, we know that if we want 
gain the connecting hose’s operating characteristics and 
its effect on the cooling performance of the cooler, the 
main target parameter, p H , must be determined ac¬ 
cording the dynamic performance simulation of the 
cooler operation, which was described in detail in Ref. 
[3]. p H is one important boundary parameter for gas 
pressure along the connecting hose to be determined. As 
for the inner nodes, their parameters can be calculated 
by RFEP and VM (Reverse Forward Euler Prediction 
and Verification Method) [16]. 

Given the cooler dimensions, the material physical 
properties, the operating conditions and the working 
fluid properties, we carry out the calculation on the 
pressure along the connecting hose according to the flow 
chart in Fig. 3. In order to boost computational accu¬ 
racy and convergence rate, the mass flow rate through 
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Adjust given conditions 



Fig. 3. Flow chart for connecting hose characteristics and cooler performance analysis. 


the convergent nozzle is verified using the following ^ ( )/ c + AF C sin( 1 . 847 t/t) {x c ^ 0) 

formula with a correction factor Vc = {v c + A Vc sin(2.1 6nfx) (x c < 0) 


(30) 


m = /.iA t 


2 7 P f Pp 

y + 1 V V p 


2/r 


( Pi 

\p 


( 7 + 1 )// 


(29) 


/< is a correction factor, which is determined by the ex¬ 
perimental data. 

The measured volume variations versus time for the 
compression and cold chambers deviate the assumption 
of the sinusoidal volume variations for the compression 
and cold chambers. So we utilized equations for the 
fitted curve to experimental results in place of Eq. (20) 
and (22), respectively. According our experimental re¬ 
sults [2,3] while coldtip temperature maintaining 80 K, 
the rectified displacement equations can be written by 


f V H - AVk sin(1.847t/T + 0.92<2>) (i D ^ 0) 

\F h - AE H sin(2.167c/T+ 1.08<2>) (jc d ^ 0) 

(31) 


The cooler performance simulation is adopted the 
method described in the Ref. [3], The iteration results 
should be satisfied by periodic stable conditions for all 
time-dependant parameters and the energy balance 
condition for the matrix of the regenerator during a 
cycle. The convergence criteria is: 


|fl(T)-0(t + r)| 

m\ 


( 32 ) 
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5. Predicted results and comparison with experiment 

Verification of our numerical modeling was carried 
out by comparing the simulation results to experimental 
data. The tested machine is a one-stage Oxford split- 
Stirling cryocooler constructed at Cryogenic Laboratory 
in Shanghai Institute of Technical Physics (SITP). The 
main parameters are as follows: compressor displace¬ 
ment 1.6 cm 3 ; compression piston stroke 8 mm; displacer 
stroke 4 mm; frequency 40 Hz; initial charge pressure 1.3 
MPa. The ambient temperature is 303 K and the coldtip 
temperature is 80 K. Dynamic pressures in the com¬ 
pression and hot chambers were measured with two pi¬ 
ezoelectric transducers and charge amplifiers. The 
displacements of the compression piston and the dis¬ 
placer were measured with two linear variable-differen¬ 
tial position transformers, enabling the volume change to 
be determined according to the driven mechanism re¬ 
spectively. The coldtip temperature was measured with a 
platinum resistance thermometer. All experimental data 
were reduced and compared by computer program. All 
calculated and tested hoses made of cooper and with a 
wall in thickness of 0.5 nun, are manufactured by Ji- 
angyin Copper Tube Corporation. 

Figs. 4 and 5 show connecting hose inner diameter vs. 
cooler performance, and connecting hose inner diameter 
vs. pressure phase lag and amplitude attenuation for a 
given hose length of 220 mm, respectively. Note the 
phase lag and amplitude attenuation are in comparison 
between the overall simulation and experimental pres¬ 
sure waves in this paper. The experimental results are in 
good agreement with corresponding ones on the calcu¬ 
lated data fitted curves in these two figures. This fact 
indicates our simulated program is feasible and in good 
accuracy. In Fig. 4, with the connecting tube inner di¬ 
ameter decreasing, the cooler net cooling power reduces 
and the compressor input power increases sharply, re¬ 
spectively. The pressure wave phase lag and amplitude 



Fig. 4. Connecting hose inner diameter vs. cooler performance. 



Fig. 5. Connecting hose inner diameter vs. pressure phase lag and 
amplitude. 

attenuation effects caused by the connecting hose were 
clearly reflected in Figs. 5 and 8. Although the pressure 
fluctuations in different chambers for a cooler with a 
smaller inner diameter of connecting hose are increased, 
to some extent, than the corresponding ones for a cooler 
with a larger inner diameter of connecting hose; the 
phase lag and amplitude attenuation for the pressure 
waves increasing sharply with the decreasing inner di¬ 
ameter of the connecting hose, were one cause for the 
cooler net cooling power reducing and compressor input 
power increasing sharply. In the meantime, we combine 
variation of mass flow rate and dynamic pressure in cold 
chamber versus time for cooler motion phase difference 

= 78 ° and coldtip temperature 80 K, which were given 
in Fig. 9. Although the mass flow rate wave and the 
pressure wave for the gas in the cold chamber have larger 
amplitudes for the cooler with a connecting hose (length 
220 nun, inner diameter 1.5 mm) than those for the 
cooler with a connecting hose (length 220 mm, inner 
diameter 2.0 mm), the phase difference between the mass 
flow rate wave and the pressure wave for the gas in the 
cold chamber are increasing for the cooler with a con¬ 
necting hose of a smaller inner diameter. Furthermore, 
the increasing phase difference between the mass flow 
rate wave and the pressure wave for the gas in the cold 
chamber directly contributes more in lowering the cooler 
cooling power output (i.e. cooler net cooling power). 

Figs. 6 and 7 show connecting hose length vs. cooler 
performance, and connecting hose length vs. pressure 
phase lag and amplitude for a given hose inner diameter 
2.0 mm, respectively. With the connecting hose length 
increasing, the cooler net cooling power and the com¬ 
pressor input power reducing, but the former has more 
large slowdown rate than the latter. So it is with the 
phase lag and the pressure amplitude attenuation. 

Fig. 8 shows pressure variations versus time for the 
three chambers while cooler using 220 mm long 
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Fig. 6. Connecting hose length vs. cooler performance. 



Con. Hose Length / mm 


Fig. 7. Connecting hose length vs. pressure phase lag and amplitude. 



connecting hose of inner diameter 2.0 mm (Fig. 8A) and 
1.5 mm (Fig. 8B) respectively. The simulated dynamic 
pressures for gas in compression chamber, which were 
derived from Eq. (20), (22) and Eq. (30), (31) respec¬ 
tively, show that simulated result by utilizing rectihed 
displacement Eq. (30) and (31) has higher accuracy to 
the corresponding experimental result. In Fig. 8A, The 
tested pressure amplitudes for the compression and hot 
chambers decrease by 11.51% and 14.8% of their simu¬ 
lating results respectively; the tested pressure amplitude 
for hot chamber decrease by 12.0% of the tested result 
for the compression chamber. Comparing with the 
simulating result for the compression chamber the simu¬ 
lating results for the hot and cold chambers lag 0.5 ms 
(7.2°) and 0.78 ms (11.3°) in phase respectively. The 
tested results are in phase with the simulating results for 
the compression and hot chambers. The tested maxi¬ 
mum pressures for the compression and hot chambers 
are in consistent with the simulating ones, however, the 
measured minimum pressures increasing by about 4.8% 
and 6.0% than their simulating results respectively. All 
these pressure amplitudes attenuation and phase delay 
occur analogously in Fig. 8B for the cooler using a 
connecting hose of 220 mm in length and inner diameter 
1.5 mm, however, the corresponding values are reduced 
because of the increased compressor pressure fluctuation 
(the simulated (measured) pressure amplitudes are 307 
(283) kPa and 277 (247) kPa for cooler with 220 mm 
long connecting hose of inner diameter 1.5 and 2.0 mm, 
respectively), and pressure ratio (the simulated (mea¬ 
sured) pressure ratios are 1.60 (1.54) and 1.46 (1.38) for 
cooler with 220 mm long connecting hose of inner di¬ 
ameter 1.5 and 2.0 mm, respectively). Besides all above, 
Fig. 8A and B indicate that time of compressing process 
for cooler with equally long connecting hose of smaller 
inner diameter is longer in a working cycle. This phe¬ 
nomenon was caused by more load force and compres- 



Fig. 8. Variation of pressure with time at coldtip cold temperature 80 K and 0 = 78°, connecting hose length 220 mm. 
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sor piston providing with more work on gas for the 
cooler with a connecting hose of smaller inner diameter. 
By comparing Figs. 8 and 9, we can deduce that with the 
connecting hose inner diameter lowering down, the 
phase differential angle between the mass flow rate and 
pressure variations for gas in the cold chamber increases 
from about 10.2 ms (147°) for the cooler with connecting 
hose of 220 nun in length and 1.5 mm diameter to 11.5 
ms (165°) or so for the cooler with connecting hose of 
220 mm in length and 2.0 mm diameter; this increasing 
phase differential angle was the main cause which gave 
rise to the net cooling power rapidly lowering down. The 
above phase lag and amplitude reduction reflected in 
detail the damp effect of the connecting hose. 

Fig. 9 shows variation of mass flow rate and dynamic 
pressure in cold chamber versus time for the cooler with 
different inner diameter connecting hose. Connecting 
hose length vs. cooldown time from ambient tempera- 



Fig. 9. Variation of mass flow rate and dynamic pressure in cold 
chamber versus time. <t> = 78°, coldtip temperature 80 K. 1, 3: con¬ 
necting hose inner diameter 2 mm; 2, 4: connecting hose inner diameter 
1.5 mm. 
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Fig. 10. Connecting hose length vs. cooldown time from ambient 
temperature to coldtip 80 K for the connecting hose of inner diameter 
2 mm. 



ture to coldtip 80 K for the connecting hose of inner 
diameter 2.0 mm is given in Fig. 10. The experimental 
results of cooldown time are approximately in propor¬ 
tional with the connecting hose length. 


6. Conclusion 

Combined with dynamic modeling for a one-stage 
Oxford split-Stirling cryocooler, connecting hose’s op¬ 
erating characteristics and its effect on the cooling 
performance of an 80 K Oxford split-Stirling-cycle 
cryocooler has been studied systematically in this paper. 
By our comparing the simulation and experimental re¬ 
sults, the connecting hose’s influence of cooling power 
reducing and compressor input power increasing for the 
Oxford split-Stirling-cycle cryocooler is due to the phase 
lag and amplitude attenuation effects of the pressure 
wave propagation through the connecting. The inner 
diameter of the connecting hose should not less than 2 
mm because it has relatively insensitive influence to the 
performance of the cryocooler in this case; the length of 
connecting hose should be as short as possible under the 
practical constraints because of its lowering the cooler 
performance and prolonging the cooldown time. Ac¬ 
ceptable agreement between experiment and simulation 
has been found. The developed computer program can 
be utilized to select the suitable connecting hose di¬ 
mensions, furthermore, develop an understanding of the 
internal process, and predicting the dynamic character¬ 
istics and operating performances. 
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